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SUM4ARY 


A method of calculating the temperature of thick walls has been 
developed in which are used relatively new concepts, such as the t ime 
series and the response to a unit triangle variation of surface temper- 
ature, together with essentially standard formulas for transient tem- 
perature and heat flow into thick walls . The method can be used without 
knowledge of the mathematical tools of its development. The method is 
particularly suitable for determining the wall temperature in one- 
dimensional thermal problems in aeronautics where there is a continuous 
variation of the heat-transfer coefficient and adiabatic -wall tempera- 
ture. The method also offers a convenient means for solving the inverse 
problem of determining the heat-flow history when temperature history is 
known. 

A series of diversified problems were solved by exact analysis as 
well as by the new method. A comparison of the results shows the new 
method to be accurate. The labor involved is very modest in considera- 
tion of the nature of the thick -wall temperature problem. T.lml+.-tn g 
solutions for the "infinitely thick" wall and for walls so thin that 
thermal lag can be neglected were also obtained. 


INTRODUCTION 


In aeronautical applications, external surfaces are heated by the' 
impact and friction of the air. For cases in which the structural tem- 
peratures never reach equilibrium, the transient temperatures of the 
surfaces often govern the design; and it is necessary to be able to pre- 
dict these temperatures. 

Literature on transient temperatures in thick walls dates from the 
classical works of Fourier. Perhaps the most extensive work on the 
subject is given in reference 1. Most literature giving the solution to 
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the transient temperatures in thick walls is based on the premise that 
the temperature history of one or more principal surfaces is known or 
given. Only a l imi ted amount of literature is available relative to 
transient temperatures in thick walls under the influence of forced 
convection. The forced-convection equation for heat transfer in aero- 
nautical applications is q = h(T aw - T), which states that the rate of 
heating q is proportional to the difference between the adiabatic -wall 
temperature T aw and the wall temperature T. The coefficient of pro- 
portionality is the heat-transfer coefficient h. In the classical 
problem of the convection heating of a thick wall, h has been assumed 
to be constant. In the usual aeronautical application, the fact that h 
varies with time is the source of the difficulty in obtaining a solution. 

The thick -wall case treated in this paper is the one governed by 
Fourier's classical partial differential equation 


dT _ v 
St 


c£t 

Sx 2 


In the case governed by this equation, the wall is composed of a homo- 
geneous material and the temperature gradients and heat flow parallel to 
the surface are negligible. In one boundary relation for this case, the 
convective heat rate is equated to the heat absorbed by the wall or to 
the product of the conductivity and the temperature gradient in the wall 
at the heated surface; that is. 


- T > ■ K (S) W 


Since h occurs as a product with T in this boundary equation, the 
usual procedures of operational calculus do not apply. When solutions 
for the temperatures of thick walls have been necessary in aeronautical 
work, the method generally used has been to divide the thick wall into a 
number of slabs in order to make a step-by-step numerical Integration of 
Fourier's equation of heat flow. Since steps in both distance and time 
must be taken, the procedure is tedious and time consuming unless the 
use of a high-speed automatic computing machine is resorted to. If it 
is necessary to do the work without the use of such equipment, a method 
introduced by Schmidt (ref. 2) wherein some of the calculations are 
accomplished graphically may be used to reduce the labor to some extent. 
This method is known as the Schmidt plot method. 
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In the present paper a simple method is developed for the calcula- 
tion of the temperature history of the surfaces of a thick wall or of 
any plane within the wall. The procedure is to select from a table a 
set of coefficients which depend on the physical properties of the wall. 
These coefficients and other data are substituted into explicit algebraic 
formulas to determine the temperature of the heated wall surface. If the 
heat-transfer coefficients are known, no guess or iteration procedure is 
required. As can be seen by the results of the example problems presented, 
the accuracy can be as good as is desired. For equal time -step sizes, the 
method is more accurate than more laborious numerical methods. 

The simplicity of the results depends on two factors: One is the 

suppression of the variable x representing the distance into the wall by 
using an integrated form of Fourier's equation and assigning a value of x 
corresponding to the heated surface. The other is a mathematical device 
known as the time series introduced by Tustin (ref. 5) ■ The time series 
is defined in appendix A. Reference 3 also introduced various manipula- 
tions of the series. The multiplication of two series is an important 
manipulation by means of which specific results can be generalized. Other 
writers (ref. 4, for example) have also presented various manipulations of 
the series. 

The present paper is divided into two parts, analysis and application. 
The section on analysis includes a treatment of the determination of the 
temperature history for the special cases of the thermally thin wall and 
the infinitely thick wall as well as for the wall of intermediate thickness. 
The inverse problem of determining the heat flow corresponding to a known 
temperature history is also discussed. Although the 'method was set up for 
the purpose of predicting wall temperatures in engineering applications, it 
has also been found to be suitable for research applications wherein the 
transient skin temperature is measured and the heat-flow and heat-transfer- 
coefficient histories are deduced. Appendix A gives background material 
pertaining to the use of time series that may be an aid to a study of the 
analysis. Appendix B gives a summary of analytical temperature and heat- 
flow formulas used either as a basis of analysis or used in the solution 
of examples to test the accuracy of the present method. In the section on 
application the computing formulas are reviewed and several examples of 
their use are given. Because of the explicit nature of the temperature 
formulas, it is not necessary to study the analysis to use the results; 


SYMBOIS 


Aq, A-j_, Ag, . . . A m dimensionless coefficients to determine heated- 

surface temperature history 

slope of wall surface temperature with respect to t ime 
°F/hr J 


b 



k 
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Bo, ®i, Bg, 

c 

e 

P(x,t) 

F A (t) 

G 

Tih 

h 

H 

k 

K 

l 

^ 2 , m 3 
m 
n 

q. 

Sn 


. . B m dimensionless coefficients to determine inside- 
wall temperature history- 

specific heat, Btu/lh-°F 

base of natural logarithms 

any function corresponding to reference slope y = (l/5)t 
any function corresponding to arbitrary . control line 

y = y(t) 

any function corresponding to unit triangle control line 

heat capacity of wall, pel, Btu/(sq. ft) (°F) 
heat-transfer coefficient at_ x = 0, Btu/(hr)(sq. ft)(°F) 
heat-transfer coefficient at x = l, Btu/(hr)(sq_ ft)(°F) 
heat -transfer number, h&jt^/l6G 

heat -transfer number for 'infinitely thick wall, 3h^KkB/8K 

diffusivity, K/cp, sq. ft/hr 

conductivity, (Btu) (ft)/(hr) (sq. ft)(°F) 

wall thickness, ft 

memory coefficients, dimensionless 

term designating time in multiples of basic interval 5 
term number in infinite series 

instantaneous heat -transfer rate due to uniform tempera- 
ture rise of heated wall surface of 1° in time 5, 
Btu/(hr)(sq. ft) 

average heat- transfer rate from time (m - 1)5 to mS 
due to uniform temperature rise of heated wall surface 
of 1° in time 6, Btu/(hr)(sq. ft) 
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heat-transfer rate corresponding to unit triangle variation 
of surface temperature 


^A,m 


r 

R 


t 

T 


average heat -transfer rate from time (m - 1)5 to mS due 
to unit triangle reference temperature variation of 
heated-wall surface, Btu/(hr)(sq. ft) 

heat-flow history at heated surface due to unit temperature 
step of heated surface 

radiation rate, Btu/(hr)(sq. ft) 

radiation term, r&it^/l 6 g, °F 

radiation term for infinitely thick wall, 5R^jtk5/8K, °F 
time, hr 

heated-wall-surface temperature, °F 


Tf, Tq, . . . 0^ T expressed as time series, °F 


T]_ value of step in wall -surface temperature, °F 

T aw adiabatic -wall temperature or effective boundary-layer 

temperature, °F 

temperature of inside (unheated) surface or of any plane 
within wall, °F 

V(x,t) temperature response to unit step in T aw , °F 


x distance through wall, ft 

y ordinate of control line or altitude of triangle 

Oq positive roots of auxiliary equation in analytical 

solution of wall temperature , 

8 basic time interval in time series 

tj ratio of heat-transfer coefficient at cooler wall surface 

to heat-transfer coefficient at heated wall surface 

difference in temperature between heated surface and any 
other plane due to unit triangle variation of heated 
surface, °F 


e 
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e r difference in temperature between heated surface and any 

other plane due to uniform reference -temperature rise 
of 1° in time 5, °F 

p weight density, lb/cu ft 

r dummy time variable, hr 

A prime denotes the derivative with respect to time. 


ANALYSIS 

TEMPERATURES ON OUTSIDE SURFACE 


Problem 

The wall considered in this paper is composed of a homogeneous 
material, and the temperature gradients and heat flow parallel to the 
surface are negligible. The transient temperatures of the heated or 
outside surface of the wall are determined by means of Fourier's equa- 
tion which governs the heat flow through the wall: 


ST _ k S^T 
St ' Sx 2 


( 1 ) 


The wall properties expressed by the diffusivity k are constant. The 
boundary conditions are given by the statements that the flow of heat at 
the unheated face of the wall (where x is taken as zero) is zero, that 
is. 



( 2 ) 


and that the heat transferred to the heated face is given by the relation 


h(T aw - T) = K(dT/ax) x=l 


(3) 
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For convenience , the initial tempe rat nr e is taken as zero at zero time : 


T = 0 


(t = 0) (4) 


Since both h and T are functions of t and occur as a product in 

equation (3), the transform procedures of operational calculus do not 
apply. The problem may be stated in another form by means of an integral 
equation expressing the heat balance at the heated surface. 

Let ig(t) be the heat-flow history at the heated wall surface at 

x = l due to a unit step in that wall -surface temperature at t = 0. 

Then the heat flow q(t) due to a temperature history T(t) ' of the 
heated surface that is continuous and is zero when t = 0, but is other- 
wise an arbitrary variation, can be expressed by Duhamel's integral in 
the form indicated in equation (5). In this integral T' represents 

— T(t) and r is a dummy time variable, 
dt 

q.(t) = T*(t - T)q B (r)dT (5) 

A heat balance is formed at x = l by equating the integral in equa- 
tion (5) to the left member of the boundary condition expressed by 
equation (3) J 


M T aw ' T ) = 



T'(t - r)q s (T)dT 


(6) 


The determination of T from equation (6) establishes the solution. 


Method 

The wall temperature T is determined, in general, from equation (6) 
for thermally thin, thick, and Infinitely thick walls . The method first 
discussed is based on thick walls, and modifications of this method 
are introduced for the special cases of thermally thin and infinitely 
thick walls. 
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Thermally thick valla . - 

Time series: In the calculation of wall temperature T for thick 

walls, the right member of equation (6) is replaced by the product of 
two time series. As explained in detail in appendix A, a time series 
is the value (here the ordinate) of a function of time at successive 
equal increments of time 5. Thus, any variation of wall surface tem- 
perature can be expressed as the series 


T - T-j_, Tg, T t )) . . . 1^ 


(7) 


When a wall surface has a unit triangle variation of surface temperature, 
the surface temperature increases from 0° at a constant rate to a value 
of 1° at the time t = 8 and decreases at a constant rate to the value 
0° at t = 28. The heat flow corresponding to a unit triangle variation 
of surface temperature can be expressed by the series 


^A q A,P q A,2' q A,3' 


A A,m 


( 8 ) 


The product of equations (7) and (8) gives the instantaneous values of 
heat flow due to the temperature variation given by equation (7) and 
could be used to eliminate the integral in equation (6). However, a 
slight variation or refinement of the method is introduced which has 
been found to increase the accuracy of the results without increasing 
the labor involved. 

If m represents the term number of a time series, the product m5 
represents the corresponding time. The refinement consists in averaging 
the heat flow from the time (m - 1)8 to the time m&. Let q^ m 

represent the average heat flow over this interval due to a triangular 
variation in surface temperature. Then the heat-flow history due to the 
triangular temperature variation can be represented by the series 


q A %1' q A,2 ; q A, 5' ’ ‘ ’ q A,m 


(9) 


where 1 is the average for the time 0 to 8, q^ is the average 

for time 8 to 28, and so forth. 
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The heat-flcrw history expressed as the average heat flow over 
successive increments 5 , hut due to the temperature variation (7), is 
given by the product of equations (7) and (9): 


S - (?!_, T 2 , *5, . . . §A } 2> ‘ ‘ * ^A,m) ( 10 ) 


Such a multiplication actually gives the result by forming the proper 
superpositions, as demonstrated in appendix A. 

In order to adjust the boundary condition expressed by equation (3) 
or the left side of eq.ua tion (6) to represent an average flow of heat ^ 

over the interval 8, the average flow of heat from the boundary layer is 
approximated by the mean of the values at the beginning and end of the 
interval. Thus, for the interval ending at t = mS, has the value 

^m = 2{^^ aw ~ "^m + ^m-l(^aw “ (H) 

If radiation is important to the problem, it may be included. Let 
r m be the rate of heat radiation per unit area at the time m 5 . Equa- 
tion (ll) may then be written 

* 4 n = 2 “ -Hm r m + ^m-l(-^aw “ ^)m-l ” r m-l] (H & ) 


With or without radiation, the heating history, or variation of q, 
may be expressed by the series 


5 Qq, <Ly . • • ^ 


( 12 ) 


wherein each term has the value given by equation (ll) or (lla). Equating 
the average heat flow given by equation (12) to that given by equation (10) 
yields 


0 - 2 ’ 




-(* 


1 ' 2 : 






( 13 ) 
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In order to evaluate equation (15)* "the product in the right-hand 
member is expanded by algebraic multiplication and values of both mem- 
bers of the equation pertaining to equal time Eire equated. The following 
set of equations results (for simplicity, the radiation terms are not 
written in) : 


|[hl(T a w - T ) L + *o( T av " T ) 0 ] = " 

2 j^^aw “ 1)2 + ^1 (^aw l) l] = ^2^A,1 + 


> (14) 


\ ^m( T aw “ T ) m + ^-1 ( T aw “ T ) m _i 


= + ^-1^A,2+ • • • +T l^\,m| 


Equations (14) can be rearranged to obtain the equivalent equations for T: 


T 1 - 


ll T aw , 1 + Ip^aw , 0 
2 ^A,1 + h l 


^ __ , 2 + ll^aw , 1 h L T x - 2X^,2 

2 2 ^A,1 + h 2 


> (15) 


•*TE = 


M, * ( bT a« - M ) m -i - 2 (9a,g^-i * • • • + 3A,m T ]) 

£ 5a,i + "m 


The values of q^ m must be derived. 

Determination of average heat flow due to unit triangle variation 
of surface' temperature : In order to obtain the average heat flow due to 

a unit triangle variation of surface temperature, the average heat-flow 
rate ^ due to a uniform increase in the surface temperature of 1° in 

each unit time 6 must be obtained. The average is taken over the 
time 6 by integrating the instantaneous heat-flcw rate from t = (m - 1)6 
to t = m5 and dividing by 6. Ihis determination is carried out in 
appendix B. The result is 
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(16) 


In equation (l6), and are the summations 


■V 


-£ 


-m(2n-l) 2 £ S§ 


n=l js£ k£ (2ll . x) 4 
4 Z 2 V 


Aa-l “ 


e 4 12 


n=l it 


£=■ — (2n - l) 2 

T Z 2 V 


(IT) 


and G is the heat capacity of the wall per square foot per °F and is 
the product of weight density, specific heat, and wall thickness: 


G = pcZ 


(18) 


The average heat flow due to a unit triangle variation of surface 
temperature q^ m is obtained by the superposition of the heat flows 

q due to three linear variations of wall temperature as follows (for 
further details, see the development of equation (Al) in appendix A) : 


^A,m " ^ < 3m-l + ( 3m-2 


( 19 ) 


Expanding equation (19) by substituting for q^ from equation (l 6) gives 






(20) 
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Substituting in equation (20) successive values of m, starting with 
m =1, ignoring any terms with negative subscripts, and collecting terms 
result in the following equations : 


3a,i 

- 8G ^ 2 + 

9a, 2 

- 80 ( id 
* 

9A,3 

' J|( A 3 - 

9A,m 

. 80 /, . 
*28 ^ 

the 

quantities 


( 21 ) 


quantities in parentheses usually retain significant values after the 
completion of the temperature triangle which created them. For this 
reason, in accordance with the notation of reference 5, they are called 
memory terms and are designated by the symbol M. With this notation, 
equations (21) become 


*25 1 


^,2 = «2 




- 8G m 

%>« - n 2g *» 




(21a) 


Hence, 


the following equations for memory terms are established: 
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Ml = ^ + A L - Aq 




M 2 = - + Ag - 3A]_ + 2A 0 

» Aj - 3^2 + 3 Ai - Aq 




(22) 


Mm - Am ^Aju^ + 3A m _2 - A m _j 


Obviously, the values of M are combinations of the values of A. 
However, as is explained in greater detail in the section entitled 
"Application," it is usually not necessary to calculate the values of 
A and M since the values of M listed in table I may be used. The 
value of M decreases with increasing term number and sooner or later 
further terms can be neglected. 


Resulting temperature formulas : Equations (21a) give the value of 

q sought to complete equations ( 15 ) ; therefore, equations (21a) are 
m 

substituted into equations ( 15 ) . Hie result can be simplified by 
dividing through by l6G/rt 2 6 and letting 



If radiation is important, the appropriate terms are included by using 
equation (lla) rather than equation (ll). Since equations (15) are 
being divided through by l6G/rt 2 5, the radiation term R is defined as 


R 



(24) 


With the substitutions of equations (23) and (24) in equations ( 13 ) , 
the final results, including terms for radiation, are 
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Ti 


^l^aw,l + ^0^aw,0 ~ -^1 ~ ^0 
M X + H x 


t 2 = 


= ^V, 2 + H l T aw,l - H 1 T 1 ~ M 2 T 1 - *2 - R 1 

Mi + H2 


M 25 ) 


_ + (^aw ~ ^m-l " **2^-1 ~ ^^-2 " • • • ~ Vl ~ ~ ^-1 

^ . M L + 


Infinitely thick vails . - 

General considerations: If a vail is thermally very thick and is 

heated rapidly so that the uhheated side experiences little heating, it 
is convenient and accurate to assume that the wall is infinitely thick. 
The same formulas, equations (25 ), are used to compute the wall surface 
temperature . However, instead of the values of M for a particular 
wall or diffusion number, the values of M which are used are always a 
fixed set of numbers which are now derived. The values of H and R 
are also changed. 

Determination of average heat flow due to unit triangle variation 
of surface temperature: The determination of the heat flow due to a 

unit triangle variation of surface temperature of an infinitely thick 
wall depends upon the instantaneous heat flow into the surface due to a 
unit rise in surface temperature in unit time. From page 110 of* ref- 
erence 4 , the instantaneous heat flow is equivalent to 2 K/tT/V 5 k. Since 
the heat-flow rate is proportional to the surf ace -temperature slope, the 
instantaneous heat transfer due to unit rise of surface temperature in 
the time 8 is 2 Kv/t/s/rtk. This expression is integrated with respect 
to t between the limits (m - 1)8 and m8. Dividing by 8 gives, 
for the average heat -flow rate ^ over the interval 8 terminating 

at m&, 

= d«L.[ m 5/2 - (m - 1)3/2] 

5/SS 

The us uhl superposition required to change the result of the slope func- 
tion to that of the unit triangle input function is accoirplished by sub- 
stituting equation (26) into equation (19): 
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QA, m = ^“=jm 5/ 2 - 3 (m - l ) 5 / 2 + 3 (m - 2) 5 / 2 - (m - 5 )^ 2 ] ( 27 ) 

Resulting temperature formulas . - If the bracketed quantity in equa- 
tion ( 27 ) is designated as the memory coefficient then 


= m 5 / 2 _ 3( m _ 1)5/2 + _ 2 )^/ 2 - (m - 3)^ 2 (28) 


A dimensionless heat-transfer coefficient (suggested by eq. (27)) is 
defined as 


3hy/itk6 

ac 


(29) 


and a corresponding term for radiation is defined as 


R 


00 



(50) 


The substitution of equation ( 27 ) into equations ( 15 ) again results in 
equations ( 25 ). Hence, equations ( 25 ) are used to obtain the heated- 
surface temperature of the infinitely thick wall as well as of walls of 
intermediate thickness, except that M, H, and R for infinitely thick 
walls are defined by equations ( 28 ) to ( 30 ). Inspection of equations ( 28 ) 
to ( 30 ) indicates that the wall material properties and time-step size are 
expressed by equations ( 29 ) and ( 3 0 ), while the memory terms are invariant 
with wall properties or step size. Substituting successive integers for 
m from 1 to 20 into equation ( 28 ) gives the following corresponding 
values of M: 


II 

H 

• 

O 

M6 = 

-0.040234 

Mjjl = - 0.012874 

Mj_g = -0.006807 


Mg = -O.I71573 

“r - 

-0.029556 

m 12 = -0* 011069 

M 17 = -0.006157 


m 5 = -0.289129 

MS = 

-0.022885 

m 15 = -0.009650 

m i8 = -0.005605 

’ (51) 

M4 = -O.IO3I76 

My = 

-0.018412 

= -O.OO85H 

= -0.005130 


Mcj = -O.O5963O 

& 

O 

tl 

: -0.015232 

m 15 = -0.007580 

Mg 0 = -0.004719 
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These values of M, along with values of E^, and from equa- 

tions (29) and ( 30 ), can be used in all problems wherein the wall is so 
thick relative to the heating rates and times involved that the wall 
behaves as though it were infinitely thick. 

Thin walls . - When a wall is thermally thin, the temperature drop 
through the wall becomes negligible and the problem is simplified by 
assuming that all interior temperatures are equal to the surface tem- 
perature. The heat absorbed by the wall during any time interval 5 
must be equal to the gain in enthalpy or total heat during this time. 
Hence, 


^m ~ g^5n^m “ ^m-l-^m-l] 


(32) 


Equating the average rate of gain jr enthalpy as given by equation ( 32 ) 
to the average rate of heat transfer through the boun dar y layer as given 
by equation (lla) results in the following heat balance: 


^m^aw - T) m - r m + hjn-l^Tav - T^ m _^- r^-i 



(33) 


Solving for gives 


T X =• 


T 2 =* 


m 


h l T aw,l + ^0^aw,0 " *^0 + § G 0 T 0 " r l _ r 0 



iT 1 

+ 

< 5 * 

OJl to 


h 2 T aw,2 + 

h l T aw,l “ h l T l + | 

CJlTi - r 2 - r x 


| G 2 + h 2 


^m^aw,m + 

(hT av - hT + | GOA 

- r m - r m _i 


^ 0 'm-1 


g + ^m 


(34) 
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Any variation of G with temperature is accounted for by equa- 
tions (34) . If the wall properties do not change over the temperature 
range covered, obviously, G is a constant. If G Is considered to be 
constant, equations (34) can be derived from equations ( 25 ) as follows: 

As the diffusion number k&/z2 becomes large, all values of A approach 
the value * 2 /8, and from equations (22) it Is seen that the only memory 
terms not identically equal to zero are M]_ and Hq, which have the 

values rt 2 /8 and -it 2 /8. Eliminating M from equations (25) and uti- 
lizing definitions (23) and (24) yield equations (3*0* 

In equations (34) the terms h Q , and Tq have been 

retained since, unlike the thick-wall problem, it is here convenient for 
T 0 to have any value. These equations have considerable advantage since 

t he need for temperature extrapolation is reduced, if not eliminated. The 
equations tend to give accurate results and, as is shown subsequently, are 
suitable for the use of relatively large time Increments . 


INS HE TEMPERATURES 


If the heated wall surface is called the outside surface, the temper- 
atures at other parallel planes may be called inside temperatures. In 
particular, this paper is concerned with the inside surface temperature. 
According to the notation used in this paper, the Inside surface 1 b desig- 
nated by x/l = 0, the outside surface by x/l = 1, and other planes by 
values of x/l between 0 and 1. 


Consider a wall, initially at zero temperature, which has the heated 
surface x = Z raised at a reference temperature slope T = (l/6)t, while 
the surface x = 0 is insulated. The difference In temperature 0 r 

between the heated surface and any plane x is shown by equation (Bll) of 
appendix B to be 


e r = If. 
r *5 k5 


£ co-s 

fel (2n - 1)5 


(2n 




n=l 


(-D 11 - 1 

(2n - l) 5 


cos (2n 


“(2n - 1)| f 


-(2n-l) 2 £ ££ 
^ Z2 


(55) 


A set of terms Is defined to represent the summations in equation (35) 5 
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Bo . 1 if.fi cos f(2n - 1)1 5 ] 

° *3 W> Si (s£ - 1)3 L 2 1 } 


B, = (-I)"' 1 cos 

1 ks fcl (2n - 1)3 


|(2n - D| *) 


-(2n-l) 2 f 




«3 k & n=l (2n - l) 5 


cos (2n - 1)^- rl e 


-2(2n-l) 2 ^ is 


^ I- 1 )"' 1 , bob Ran 

*3 118 n=l (2n - 1)3 L 


111 Z ^ 


t]' 


-m(2n-l) 2 £ ^ 


(36a) 


If the temperature difference across the entire wall Is sought, then 
x/Z = 0 and equations (36a) "become 


B 0 

% 

®2 


16 

22 

00 

(-D : 

n-1 

«3 

k8 

n=l 

(2n - 

l) 3 

16 

Z 2 

00 

n 

(-1) : 

n-1 
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k& 
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rt 3 

k8 

n=l 

(2n - 

l) 3 



2k8 


-(2n-l) 2 ^ is 


-2(2n-l) 2 

e 



-m(2n-l) 2 


e 




(56b) 


In either case, by using equations (36a) or (36b) equation (35) may be 
written 


e r = Bq " Bm 


( 37 ) 
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Ey the usual superposition, in order to change the result of the 
reference slope function to that of a unit triangular input function, 
the temperature difference between the plane of the heated wall and the 
plane of consideration is 


0 = Bo - % - 2 (Bq - %-i) + Bq - %-2 (38) 


Assigning integral values to m and simplifying give the equations 
corresponding to successive values of time: 


e l = _B 1 + B 0 

8g = “Eg + 2B-j_ — Bq 

8 ^ = — B ^ + 2Bg — B^ 
0 ij. = -B]^ + 2Bj - Bg 


> 


0 


m 


_B m + ^m-l ~ ^ 81 - 2 ] 


( 39 ) 


A sufficient number of values of 0 for practical purposes are given in 
table I. Equations ( 39 ) mean that the temperature difference between the 
plane of the heated surface and any other plane due to a unit triangle 
temperature variation of the heated surface is given by the time series 
8^ = 8]_, 0g, 0j, . . . 8 m . In addition, any general temperature of the 

heated surface can be represented as T = T-j_, Tg, T^, . . . 0^. The 

temperature difference due to this general temperature variation is 
obtained by formally multiplying these two time series, and the inside 
temperature is obtained by subtracting the product from the heated-surface 
temperature T. Hence 


T i,l = *L " 9 1 T 1 

T i,2 = T 2 - (9l T 2 + 02^) 

T i,3 = T 3 " ( S 1 T 3 + e 2 T 2 + 9 3 T l) 

T i,m = " ^l^m + e 2 a m-l + • • • + e m T l) 

j 

Computation of m does not depend on prior computation of 
T. _, and so forth. 

1. *d. 


(*o) 


T. 


i,l J 
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If the temperature distribution through the wall is required in a 
form which is analytical with respect to x and can be differentiated 
or integrated analytically with respect to x, the procedures outlined 
in appendix C should be followed. 


CALCULATION OF HEAT FLOW FROM TEMPERATURE HISTORY 


Temperature History of Outside Surface Known 

If the heated-surface temperature history is known, the heat-flow 
history can be determined by substituting equations (21a) into equa- 
tion (10) as follows: 


5 - (Ti, t 2 , Tj, , . . ^)(m x , m 2 , m 3 , . . . mJ (41) 
Multiplication shows that the mth term is given by 


^m = 


8g 

ic^8 


(»A + “2 Vi + • 



(42) 


The heat -flow history can be determined readily from a given temperature 
history of the heated surface by means of equation (42) . The average 
heat flow over a small interval 8 can be assumed to give the instan- 
taneous rate at the center of the interval. 


Temperature History of Inside Surface Known 

If the temperature history of the outer surface or of a plane near 
the outer surface is known, the feasibility" of accurately determining 
the heat flow is excellent. If the temperature history of the inside 
surface for a thermally thin wall (kt/l^ large and hZ/K small) is 
known, it is also feasible to determine the history of the heat flow 
into the outer surface. However, if the wall is thermally thick, rela- 
tively small changes in temperatures at the inside surface may make it 
difficult to reconstruct the temperature history and heat flow at the 
outer surface. 

A rearrangement of equations (40) may be used to determine the 
outside-surface-temperature history from the inside -surface -temperature 
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history. 


Rearranging equations (40) gives 


T 2 

*m 


T i,l 


i - e L 


T i,2 + 

6 2 T 1 

1 - 

01 

T i,m + 

0 2 T m-l + 0 3 T m-2 + • • • + 0 m T l 

1 - 9 ! 


(h 5) 


The rate of heat flow at any time is then determined from equation (42) . 


APPLICATION 


GENERAL CONSIDERATIONS 


The section on application is devoted to demonstrating the solution 
of two types of problems: In the first type the heat-transfer coeffi- 

cient and adiabatic-wall-temperature histories are given and the wall- 
temperature solutions are obtained. In the second type the temperature 
history is known and the heating -rate history is computed. In each 
example, the problem chosen was a special case, selected so that its 
solution could be and was obtained by an exact analytical method. The 
degree of exactness of the present method is demonstrated by comparing 
each result with the solution calculated by exact theory. 

With one exception, radiation was a negligible consideration in the 
examples given. Although the method presented is well suited to accounting 
for radiation and includes terms for that purpose, radiation was neglected 
in all cases to make possible an exact analytical solution for comparison. 


Equations for Heated Surface 

If the temperature of the heated surface of the wall is required, 
equations (l8) and (25) to (25) are used. For convenience, these equa- 
tions are summarized as follows: 


k 


K_ 

cp 


_ hSm 2 - R fSrt 2 

i6g i6g 


G = pel 
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= H l T aw,l + ^O^aVjO ~ R 1 ~ ^0 
1 M x + H x 


T 2 


_ ^V,2 + ( HT aw ~ HT )l ~ ^1 ~ r 2 ~ R 1 


Mj_ + 


T 3 


_ H 5 T av,5 + ( m a¥ ~ ht )2 ~ M2 t 2 ~ M 5 T 1 ~ R 5 ~ ^2 

M x + H 5 


*4 


m H 4*aw,4 + fow - ~ «2^ ~ M 3 T 2 ~ M 4 T 1 - R 4 ~ *3 

Mi + 


_ + (^aw ~ ■ HT )m-l ~ ^2^1-1 ~ ^y^nx-2 ~ • • • ~ ~ ^m-l 

m " Mi + % 


Whether the objective is to compute wall temperature or to compute 
the heat flow from a known wail temperature, the first steps Involve the 
determination of the required memory coefficients. In order to minimize 
the labor involved, the recommended procedure is as follows: First, 

choose a tentative time interval 5 which seems appropriate to the 
particular problem. (A review of the examples presented herein will give 
an idea of a reasonable value . ) Then compute a tentative value of the 
dimensionless diffusion number kd/l^. From table I or II pick a dif- 
fusion number close to the one tentatively computed. The memory coeffi- 
cients M and inside temperature coefficients 8 given in the table for 
this diffusion number are to be used, and they do not therefore have to 
be computed. An adjustment in the value of 5 is made by multiplying it 
by the ratio of the tabular value of k b/l^ selected to the value of 
k b/l^ tentatively computed. Then compute from equations (l8) and (25.) 
the value of G and the values of H; if the radiation is important, 

R must be computed also (eq. (24)). The temperature history of the 
heated-wall surface is then found from equations (25). If the heating 
rate is being determined from a known temperature history, the values of 
H and R are not required. The procedure for this case is discussed 
in the section entitled "Example 6.” 

The temperature formulas were derived with the assumption that the 
initial wall temperature was zero in order to avoid writing T - Tq 
numerous times in the formula. The simplest way to handle most problems 
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is to subtract the amount that the initial vail temperature is above 
zero from both the vail and the adiabatic -vail temperatures. The last 
step in the problem is to add this amount to the solution. 


Equations for Inside Surface 


If the temperature of the unheated side of the vail is required, 
equations (40) are used. These equations are summarized for convenience: 


T i,l = T 1 ‘ 9 1 T 1 

T i,2 = T 2 " ( 9 1 T 2 + 9 2 T l) 

T i,3 = T 3 " ( 9 1 T 3 + 9 2 T 2 + 9 3 T l) 

•'■i,m “ ^th — ( 9 l^m + 9 2^m-l + • • • + ®m^l) 

Except for the case of the thermally thin vail, all equations vere 
derived for constant material properties. For small changes in material 
properties vlth temperature, it appears reasonable to use an average 
value of the properties for the temperature range involved. For cases 
in vhich material properties vary, it seems possible that a more accurate 
ansver might be obtained by varying the diffusion number or by varying 
the step size to keep the diffusion number constant; hovever, any con- 
sideration of such a technique is beyond the scope of this paper. 


APPLICATION OF METHOD IN SPECIFIC EXAMPLES 


The folloving illustrative examples vere calculated before table I 
vas prepared. Therefore, the values of the coefficients M and 0 vere 
computed for the particular vails and chosen time intervals 5 of the 
examples. The values of M and 0 used are all listed in table II, 
vhich may be considered as being supplementary to table I. 

Example 1 


Example 1(a) .- 

Problem: A copper vail vhich is l/2 inch thick is initially at a 

temperature of 0° F. One surface is heated by a boundary layer vhile the 
other side is insulated. The effective boundary-layer temperature T aw 

is initially 0° F but increases linearly at the rate of 1,000° F per 
second for 10 seconds . The heat-transfer coefficient remains constant at 
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h = 100 Btu/(hr)(sq ft)(°F). The conductivity K and diffusivity k 
of copper are taken as 


K = 227(Btu)(ft)/(hr)(ag. ft)(°F) 


k = 4.4l sq ft/hr 


Find the temperature history of both vail surfaces. 

Solution: The material properties are usually given in terms of 

the hour unit. However, since fast heating conditions may be more 
easily understood in terms of seconds, time is referred to in seconds 
and is converted to hours for use in the equations. For example, if 


6=1 sec 


1 

56OO 


hr 


then 

k5 _ (4.4l)(24) 2 _ 0 6 

TI 3600 75 


By using this dime nsionless diffusion number, the values of M in 
column 2 of the following table are obtained from table II. The values 
of T aw are listed in column 4. The value of G is given by the 

equation 

G = ¥• = / 1 wf 7 , n ' = 2.1447 Btu/(sq ft) (°F) 
k (24)(4.4l) 

The value of h =_100 Btu/(hr)(sq ft)(°F) is converted to H = 0.00800. 
The use of columns 2 and 4 in equations (25) gives the heated-wall tem- 
perature in column 5. Using the inside -surface -temperature formulas, 
equations (40), and the values of 9 in column 6 gives the values of 
Ti in column 7* 
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1 

2 

3 

4 

5 

6 

7 

Term 


Time, 

-'■aw* 


0, 

Ti, 

number 

M 

sec 

Op 

Qjt 

oj> 

op 

n 


o 

0 

0 


0 

i 

0.75160 

1 

1,000 

11 

0.580382 

4 

2 

-.35256 

2 

2,000 

36 

-.474632 

21 

3 

-.33055 

3 

3,000 

74 

-.087207 

50 

4 

-.05648 

4 

4,000 

124 

-.015291 

91 

5 

-.00990 

5 

5,000 

186 

-.002681 

144 

6 

-.00174 

6 

6,000 

260 

-.000470 

210 

7 

-.00030 

7 

7,000 

346 

-.000082 

288 

8 

-.00005 

8 

8,000 

443 

-.000015 

377 

9 

-.00001 

9 

9,000 

552 

-.000002 

478 

10 

-.00000 

10 

10,000 

673 

-.000000 

591 


The wall -surface -tender a ture curves of T and are shown in 

figure l(a). For comparison, the results calculated by the theoretically 
exact formula (eq. (BIO) of appendix B) are shewn. This formula is 


T = T. 


aw 


(44) 


where ofo represents the positive roots of the auxiliary equation 


On tan = — • (44a) 

ft. 


Setting b = 3,600,000 °F/hr and setting x = l and x = 0 in equa- 
tion (44) result in the values for outside and inside temperatures 
plotted as circles and squares in figure 1. The comparison shows that 
accurate results are obtained by the present method, the maxim um dif- 
ference between methods being about 1°. 

Example l(b) . - 

Problem: The conditions for example l(b) are the same as those for 

example l(a), except that the copper wall is 3 inches thick, or 
l = l/4 foot. 




26 


NACA IN 4105 


Solution: Since this wall is so thick, more highly transient condi- 

tions prevail throughout the heating period. A value of 5 of l/2 sec- 
ond or 1/7200 hour was therefore used. With the same procedure as used 
before, the results of the present method are given as continuous lines 
in figure l(b) , while the results from equation (44) are given by the 
symbols. The agreement is evident. 


Example 2 

Problem: Example 2 illustrates the principal advantage of the 

present method; that is, its capability of handling arbitrary variations 
of T aw and h. Neither is it necessary to know a mathematical formula- 
tion for these variations. 

A l/2-inch copper wall which is initially at a temperature of zero 
is aerodynamically heated on one side and insulated on the other. The 
time histories of T aw and h are given at l/2-second intervals in the 

following table: 


Time, 

sec 

T 

aw* 

Oj. 

, h , 

Btu/(hr)(sq. ft)(°F) 

Time, 

sec 

T 

^aw* 

Op 

, h > 

Btu/(hr)(sq. ft)(°F) 

0 

0 

36 

5.5 

5,391 

68.4 

• 5 

1,365 

41.4 

6.0 

5,356 

66.6 

1.0 

2,485 

45.0 

6.5 

5,255 

63.9 

1.5 

3,3f8 

48.6 

7.0 

5,107 

60.0 

2.0 

4,op4 

52.2 

7-5 

4,851 

55-8 

2.5 

4,620 

55-8 

8.0 

4,535 

52.2 

3.0 

4,932 

60.0 

8.5 

3,654 

48.6 

3-5 

5,119 

63.9 

9.0 

2,769 

45.0 

4.0 

5,265 

66.6 

9-5 

1,658 

41.4 

4.5 

5,345 

68.4 

10.0 

297-5 

56.0 

5.0 

5,387 

69.0 





Find the temperature history of both wall surfaces . 

Solution: The heating conditions are severe and continuously tran- 

sient, with the boundary-layer temperature rising and falling over 5 , 000 ° 
in 10 seconds . A- computing interval 5 smaller than that used in the 
first problem is therefore used. Let 5 = l/2 second = I/ 72 OO hour. The 
resulting wall-surface-temperature curves are drawn in figure 2. The 
circles and squares give the surface temperatures computed by a theoreti- 
cally exact procedure. Comparison shows the present method to be accurate. 
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Example 3 

Problem: Example 3 is the same as example 2 except that the wall is 

3 inches thick and the effective boundary-layer temperatures are slightly 
different. The initial value of T a v is 0. The subsequent values of 
T aw are given at 0.5-second intervals by the following time series: 

T aw = 1 , 365 , 2,1*84, 3,386 , 4,088 , 4,609 , 4,915, 5,094, 5 , 227 , 5 , 296 , 

5,325, 5,315, 5 , 265 , 5,1^9, 4 , 986 , 4 , 694 , 4,184, 3,489, 2 , 591 , 1 , 471 , 100 . 

Find the temperature history of both wall surfaces. 

Solution (a) - thick-wall solution: The value of 5 was taken as 

l/2 second. By following standard procedure, the wall -surface - 
temperature curves shown in figure 3 were obtained. Again the symbols 
give the results of exact theory. 

It should be noted that the "exact solution" for examples 2 and 3 
is not actually an alternate method of solution for any practical prob- 
lem but gives a solution to the particular problems only. The "solution" 
was obtained by working in reverse; that is, a heat flow was assumed and 
the corresponding boundary-layer characteristics were computed. A truly 
alternate method of solution is now considered, however. 

Solution (b) - infinitely thick -wall solution: Since the thermal 

lag of a 3 -inch copper wall is so great when subjected to the rapid 
heating specified by this problem, it appears reasonable to obtain the 
heated-surface temperature by assuming that the wall is infinitely thick. 
The memory coefficients are the same for all infinitely thick walls and 
are given by equation ( 31 ) . The same temperature formulas, the same 
values of T aw , and the same values of h are used as before, but the 

values of are given by equation ( 29 ) . The results calculated by 

this method are listed along with those from the thick-wall solution (a). 


Time, 

sec 

T, °F, calculated by - 

Thick -wall 
solution 

Infinitely 

thick-wall 

solution 

1 

12.8 

12.8 

2 

35.3 

35-3 

3 

61.3 

61.3 

4 

87 .O 

87 .O 

5 

109.1 

109.1 

6 

124.8 

124.8 

7 

132.3 

132.3 

8 

130.6 

130.6 

9 

119.5 

119.4 

10 

100.0 

99-8 
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The consistency of the alternate methods for a thermally thick wall is 
evident. The reason for the close agreement may he found in figure 3> 
which shows that the unheated surface of the 3-inch wall rose to only 
14° F. 


Example 4 

Problem: The most severe test of the present method would occur 

if there were a large instantaneous increase of T aw . While this condi- 
tion could hardly happen in flight, it might happen if a research model 
were suddenly immersed in a high-stagnation-temperature jet. Let a 
l/2-inch copper wall, initially at a temperature of zero, he instantly 
subjected to an effective boundary -layer temperature of 5>000° I* on one 
surface while no heat transfer occurs on the other surface. The heat- 
transfer coefficient is 100 Btu/(hr)(sq ft)(°F). Solve for the tempera- 
ture history of both wall surfaces for 10 seconds. 

Solution: In this case not only is there a very high tranBient- 

temperature condition initially but the instantaneous increase in T aW 

does not lend itself to approximation by the unit triangle. The 
simplest procedure is to take small steps for the first few seconds to 
minimize the errors introduced. In order to help circumvent the diffi- 
culty of calculation, an excellent method of approximating the wall 
surface temperature for the first or first few small steps is to use 
the following formula from page 109 of reference 4, which gives the tem- 
perature on the surface of an infinitely thick wall ‘for a constant flow 
of heat at the surface : 


_ ^k^aw fkt 

K V * 


(^ 5 ) 


The values of 5 used were 5 = l/8 second for 2 seconds, then 
5 = l/2 second for 8 seconds. Since the use of an equation based on an 
infinite wall is permissible for a l/2-inch copper wall for at least 
l/4 second, the values of T were computed by equation (’45) for the 
first two l/8-second steps, then by the usual equations. The results are 
presented in figure 4. Since the inside temperature T-j_ depends on the 
outside temperature T and not directly on -'■aw> there is no particular 

difficulty of approx ima tion in obtaining Tj_. Accordingly, in obtaining 
T-l, l/8-second steps were taken for 1 second to define the highly transient 
part of the curve, then l/2-second steps for the remaining 9 seconds. 

A theoretically exact solution to this problem was obtained by equa- 
tion (B7) of appendix B. The results of applying this equation are given 
by the symbols in figure 4, which shows that agreement was obtained. 
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Example 5 

Problem: A l/l 6 -inch Inconel wall is heated, by a high -temperature 

jet. If T aw = 5 , 000 ° F and. h = 50 Btu/(hr)(sq ft)(°F), determine the 
skin- temperature history for 15 seconds . Neglect radiation. 

Solution: The heat capacity of a l/l 6 -ineh Inconel wall was assumed, 

to be G = O.5229 Btu/(sq ft)(°F). The example was worked three times 
with values of 8 of 1 , 2 , and 5 seconds to show how sensitive the "thin 
wall" formula is to step size. Substituting the given constants into 
equation ( 54 ) yielded the results given by the symbols in figure 5 * For 
this example the exact theory is shown by the solid line. Evidently, 
coarse steps are permissible with this formula. 


Example 6 


Example 6(a) . - 

Problem: The temperature history of the heated surface of a l/ 2 -inch 

copper wall initially at zero temperature is given by the following time 
series in which the temperatures are separated by 1/2 second: T - 4 . 8 , 

15.2, 25.2, 40 . 2 , 58.2, 77.9, 99-4, 122.6, 146 . 5 , 169.8, 193 - 0 , 214 . 8 , 
235 - 3 , 255.2, 268.8, 281.5, 290.8, 296 . 7 , 298.0, 297 - 5 - The inside sur- 
face is insulated. Determine the history of heat flow into the heated 
surface from the given surface -temperature history. 

Solution: If the time interval used is sufficiently small, the 

average rate of heat flow over the interval is a good approximation to 
the rate of heat flow at the center of the interval. Equation ( 42 ), which 
gives the average rate of heat flow over the interval ending at t = mS, 
may be used. For a value of 8 of l /2 second, the values of M are 
given in table H. Substituting in equation ( 42 ) gives the rate of heat 
flow plotted as circles in figure 6. The solid curve gives the theoreti- 
cally exact instantaneous rate of heat flow for comparison. The results 
from equation ( 42 ) are seen to be precise. The system yielding instan- 
taneous heat flow, mentioned previously, would seem to be a natural one 
for the present problem; however, the results obtained by that system 
were found to be inferior to those presented. 

Example 6(b) . - 

Problem: The corresponding inside -surface -temperature history of 

the same wall is given by the following time series in which the tempera- 
tures are separated by l/2 second: Ti = 0.8, 4 . 4 , 11 . 7 , 22 . 4 , 56. 0, 52.6, 

71.6, 92 . 5 , 114 . 2 , 137.4, 160.6, 185.7, 205.4, 226.4, 245 . 3 , 261.6, 275.2, 

285.6, 292 . 4 , 296.5. Determine the history of heat flow into the outer 
surface by using only the given inside temperatures. 
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Solution: Equations (45) may "be used to determine the outside- 

surf ace -temperature history from the inside -temperature history. Then, 
the rate of heat flow at any time is determined as in the solution for 
example 6 (a). The factor 1/(1 - Sq) in equations (45) may be thought of 
as a magnification factor. Large values of this factor tend to cause an 
instability in the computed temperatures. In this example, if 5 is 
taken as l/2 second or 1/7200 hour, the value of kS/2 2 i S O. 5528 , 9^ 

is O. 80478 , and the magnification factor is 5*1 • On substitution in 
equations (45 ), an oscillation of period 28 builds up in the values of 
T, very slowly at first and to either side of the correct answer as a 
mean value, but in a divergent manner so that by 4 seconds the amplitude 
is 6 °. For a value of 6 of 1 second, the value of kS/Z 2 is O. 7 O 56 , 

0 p is O. 58 O 58 , and the magnification factor is 2.4. In this instance, 
an oscillation in T of period 28 and maximum amplitude of I. 5 0 
occurred. Substituting these values of T, without fairing, into equa- 
tion (42) gives the results shown by the square symbols in figure 6. If 
a larger value of 8 were used, the oscillation would be damped out but 
the accuracy would suffer because of a lack of definition of the rapidly 
varying heating rate. The particular case demonstrated is therefore a 
marginal one for the determination of heat flow from the temperatures of 
the inside surface. The instability is found to disappear for thermally 
thinner wadis and conversely to increase rapidly for thicker walls. 


CONCLUDING REMARKS 


Formulas to facilitate the determination of the transient surface 
tempera tures of thick walls from an arbitrary variation of adiabatic - 
wall temperature and heat-transfer coefficient have been developed. 
Formulas to facilitate the determination of heat flow from an arbitrary 
variation of wall surface temperature were also obtained. The numerical 
applications given demonstrate a high degree of accuracy for the present 
method. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., June 18, 1957- 
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APPENDIX A 

REVIEW OF SUPERPOSITION AND UME SERIES 


TIME SERIES 


A time series may be defined as a series of numbers or quantities 
which represent the values of a function of time at successive equal 
intervals of time. According to the notation of Tustin (ref. 3)> each 
quantity j.s separated from the others by a comma since the values corre- 
sponding to different times are not added together. The quantity at zero 
time is zero. The first quantity recorded represents the value of the 
function at the end of the first time interval. The time interval used 
is arbitrary and its size is selected to obtain the accuracy required in 
the solution of a specific problem. The symbol for the time interval is 
8. Thus, the series y = y-^, y^ s 7y • • • y m represents the values of 

the f met ion y at the times 8, 28, 36 > • • • niS. 


THE UNIT TRIANGLE 


A unit triangle is an isosceles triangle which has an altitude of 
unit magnitude and a base of 28, or two time intervals. Since 8 is an 
arbitrary time interval, the unit triangle is accordingly arbitrary. A 
plot of a unit triangle centered at t = 8 is shown in sketch 1, where 
y represents magnitude or altitude. 



Sketch 1. 



32 


NACA IN 4105 


The slopes of the sides of the triangle depend on the value of 5 and 
are equal to ±l/&. Three lines may be superimposed to represent the unit 
triangle. The equations of the three lines which may be added to repre- 
sent the unit triangle are y = (l/5)t, y = -(2/6) (t - 5), and 
y = (l/6)(t - 25). Although the triangle terminates at t = 25, the 
values of t in the equations for y can go on to infinity since the 
values of y add to 0 beyond t = 25. 


The Function Corresponding to a Reference Line 

Let F(x,t) be the solution to a boundary -value problem specified 
by a linear partial differential equation and the linear boundary condi- 
tion y = (l/5) t, where y is the value of F or one of its derivatives 
or integrals at some fixed value of x. Because of the linearity of the 
problem, the magnitude of the solution is directly proportional to the 
magnitude of the slope l/6 of -the boundary condition. For example, if 
y = (2/6)t, the corresponding function representing the solution is 
2F(x,t). Again, if y = -(2/6)t, the corresponding function is -2F(x,t). 
The slope of the line y = (l/6)t can thus be used as a reference magni- 
tude. This slope is the same as that of the left side of the unit 
triangle . 

The value of t in F(x,t) is always identical with the value of 
t in the boundary condition y = (l/6)t. Thus, if the origin is shifted 
so that y = (l/6)(t - 26), then the corresponding function is F(x,t-26). 
Particular solutions of a linear differential equation can always be 
added in linear combinations to satisfy more general boundary conditions. 
If a and b are constants, and aF(x,t) corresponds to the boundary 
condition y^ = (a/6)t, and bF(x,t-6) corresponds to yg = (b/6)(t - 6), 

then the function corresponding to the sum of the two lines y = y^ + y 2 

is F = aF(x,t) + bF(x,t-6). Let an additional property of F, as well as 
of y, be that it assumes the value 0 for any time less than 0. The 
range of time of interest is therefore from 0 to ». 


The Function Corresponding to a Triangle 

Wow - consider a function F^(t) dependent on the lines of the unit 

triangle In sketch 1 for its value just as F(t) is related to the line 
y = ( 1 / 5 ) t. In consideration of the three lines by which the unit triangle 
may be replaced, y = (l/6)t, y = -( 2 / 5 ) (t - 5), and y = (l/5)(t - 25), 
the three corresponding solutions or functions of time are F(t), 

-2F(t - 5), and F(t - 25). Because of the additive nature of solutions, 
the solution corresponding to the complete triangle may be defined as the 
sum of the solutions for the lines which compose it: 
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F A (t) = F(t) - 2F(t - 8) + F(t - 26) (Al) 

In equation (Al) each term has a value of zero for negative values of 
its argument. 


Representation of a General Curve by Triangles 

The curve A-B in sketch 2(a) is any arbitrary continuous function 
in the plane y,t which can be faired through its ordinate values 0, 
yij and y^. The curve may be well approximated by a series of 

chords such as a, b, and c if a suitable spacing 6 is used. 



(a) The curve A-B (b) The synthesis of A-B 

Sketch 2. 

Draw the lines d, e, f, g , and h of sketch 2(b) to form the isos- 
celes triangles with sides a and d, e and f , and g and h. Since the sum 
of two straight lines is a straight line, it is clear that, if the letters 
which designated the lines are now used to designate the equation of the 
line, then d + e = b and f + g = c. Of course, the line a is the first 
half of the first triangle as well as the first chord of the curve A-B. 

In designating the ordinates y^_, y^, y^, . . . y m , the altitudes of three 

triangles whose sides add up to the chords of the curve A-B are simulta- 
neously designated. The ordinate series y^, y y^, • . • y m , with 

spacing 8, is a time-series approximation of the curve A-B. In this 
case, as well as elsewhere in this report, each ordinate of a series is 
understood to represent the altitude of an isosceles triangle with a 
base width of 2B. 
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Arbitrary Control Functions 

Suppose that the curve A-B, or y(t), is to serve as an arbitrary 
reference or control function for a corresponding function Fy(t). It 

is desired to determine Fy(t) -in terms of functions corresponding to 
unit triangles, such as F A (t), given by equation (Al). The func- 
tion F A (t) is a function corresponding to a unit triangle centered 
at t = 8. If F A (t) is multiplied by y-^, there results the function 

corresponding to the triangle vith sides a and d of sketch 2(b). By 
superposition, the function corresponding to the three triangles centered 
at 8, 28, and 35 is 


Fy(t) = y L F A (t) + y 2 F A (t - 8) + y 3 F A (t - 28) (A2) 


Then, since the three triangles add up to the chords of the curve A-B, 
this Is the function, in ordinary algebraic form, corresponding exactly 
to the chords of the general reference curve y = y-^, y 2 , y^, . . . y m . 

In order to put equation (A2) in time-series form, let F A (t) = a 1 , a 2 , 
a^, . . . a^, with spacing 8. Substituting in (A2) and placing terms 
for the same time in columns and adding yield 


y x F A (t) 

= y^, 

yia 2 ; 

. y i a 3> 

- 

5 ) - 0 , 

y 2 a L , 

y 2 a 2> 

y^A^ - 

28) = 0 , 

0 , 

y 3 a i' 


Fy(t) = y x a L , (yqag + y 2 a 1 ), ( 7 ^ + y^ + y 5 a L ), ... (A3) 


The result shown by equation (A3) is obviously that which is obtained 
by formal algebraic multiplication of a^, a g , a^, ... by y^, y g , 

y^> • • 


. as follows : 
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F A (t) = 

a 2 j 

a 3> 

y(t) = y 1 , ' 

CVI 

y 5 . 

y l a l> 

y l a 2> 

y i a 3> 


y 2 a i. 

y 2 a 2 > 



y 3 V 


^y(t) - yis-i* (^1®^ + ^2 a i)^ (yi a 3 + y2®2 + y3 a i) > • • • (a4) 

Therefore , 

Fy(t) s = (j A (t) a ] [y(t) B ] (A5) 


where the subscript s denotes time-series form. 

Equation (A5) states a simple theorem which was first given in 
reference 1. It was developed with y(t) having a value of 0 at 
t = 0 and F A (t) also having a value of 0 at t = 0. If either or 

both of the series had values at t = 0, multiplication as in equa- 
tion (A5) would not be sufficient to obtain F v . 
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APPENDIX B 


SUMMARY OF ANALYTICAL TEMPERATURE FORMULAS FOR THICK WALLS 


CONSTANT FLUID TEMPERATURES 


Consider an infinite wall of thickness Z and an initial tempera- 
ture of zero. (See sketch 3*) Let the wall be suddenly contacted at 
the face x = Z by a fluid of temperature T a w while the face at x » 0 



remains exposed to a fluid of zero temperature. Let the heat-transfer 
coefficient at the face x = Z have the value h and the heat-transfer 
coefficient at the face x = 0 have the value T|h. Let the wall have 
uniform physical properties which are invariant with time. 

The flow of heat within the wall is governed by Fourier's equation 
for transient heat flow, which states that the rate of increase of 
temperature Is proportional to the rate of change of temperature gradient: 


6T = 

^ = Sx 2 


(Bl) 
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The constant of proportionality, called the diffusion coefficient, is 
equal to the ratio of the conductivity of the material to the heat 
capacity as represented by the product of specific heat and weight density 
as follows : 


k 



The boundary equation for the wall at x = l is obtained by equating the 
heat transfer between the fluid at temperature T aw and the wall at tem- 
perature T to the rate of heat transfer in the wall at x = l : 


h 




(B2) 


A similar boundary equation is written for the wall x = 0. 


Tl hT x=0 


- <S) 


x=0 


(B3) 


The initial condition is a statement that the initial temperature of the 
wall is zero: 


( T >t=0 


= 0 


(B4) 


The simultaneous solution of partial differential equations (Bl) to 
(B4) is an infinite series. 


‘ ^EEW 


I + *1 
i + n + n 


hi 


-2T< 


; av ZH 
rt=l 



eOB (f + 5^ X - BiI1 (f 4 

e 

fK 2 hi 

LET - 11 T ■ 

■ 3(1 + il) cos ctn + 


+ ') 


sin djj 


(B5) 


The parameters ct^ are angles which are the positive roots of the 
equation 
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tan an 


<%(! + q) 

a/ JL . ."hi 

^ hZ * K 


(B5a) 


Two dimensionless numbers of physical significance are present in 
equation (B5): the diffusion number kt/Z^, and the conductance or 

Nusselt number hZ/K. The first term of equation (B5) gives the steady- 
state solution or equilibrium condition. 


Special Case q = 0 

The special case of a plate that is heated by convection on one 
face while the heat transfer on the other face is so small as to be 
negligible is important. Setting the heat -transfer coefficient at the 
unheated face equal to zero corresponds to making the assumption that 
the plate is perfectly insulated at that face. If the plate is per- 
fectly insulated at the face x = 0, substitute q = 0 in equations (B5) 
and (B5 a ) and utilize the following relation to eliminate sin 


an'tana I1 = (B6) 


where represents the positive roots of the equation. These substi- 

tutions yield the following expression for temperature: 


T = T _ - 2T. 


aw 


'aw 


Cos (f °n) 

A v% 1 ir 

n=1 (- 


J2U 2 

j2 n 


1+ K- + 


JL 

hZ 


Op 2 ') cos 


(B7) 


Equations (b 6) and (BT) may be found on page 100 of reference 1. 

Special Case q = 0 and h = 00 

The special case in which q = 0 and h = a> corresponds physi- 
cally to perfect insulation on one face and a known initial temperature 
on the other face. This case is developed by setting h = 00 in equa- 
tions (B6) and (B7). Equation (B6) becomes 
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tan. = 00 


tan Oq = co 


or, since assumes the sequence 2 L, 32., 5 %., and so forth. 


on = ( 2 n - 1 )| 


Equation (B 7 ) reduces to 


T = 






n=l 


2 n 


cos 


h - «f ?]' 


-(2n-l)2 & 


(E8) 


Equation (B 8 ) may be found in standard references (ref. 4, page 19 6 , 
problem 6 , for example). In this equation, T-j_ is the step in wall- 

surface temperature used in place of T aw . 


VARIABLE FLUID TEMPERATURE 


Arbitrary Variation of T av 

Equation (B7) is an exact solution of a thick -wall boundary -layer 
heating problem which is suitable for checking the accuracy of the 
present method of computing wall temperatures . Another method, which 
is more general, can be obtained by letting the adiabatic -wall tempera- 
ture vary in a known manner, T av = F(t). Let V(x,t) be the variation 

of wall temperature due .to a 1° step in adiabatic -wall temperature . If 
the initial wall temperature is zero, the wall temperature is given by 
Duhamel's formula as 


T(x,t) 



F' (T)V(x,t-T)dr 


(B9) 


where F' is the derivative of F with respect to t and r is a 
dummy time variable. 
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Linear Variation of T aw 


In the linear variation of T aw> T aw = bt, or, in the notation of 

Duhamel's formula, F = bt and F' = b. Substituting T from equa- 
tion (B7), with T av = 1, into equation (B9) and performing the integra- 
tion give ■ 


T(x,t) = bt 


2bl2 


hi 41 -.-**) 

^ + %jr on 2 ) cos % 


(BIO) 


Equation (b 6) is used to obtain a^. 


Linear Variation of Wall Surface Temperature 

One of the fundamental equations on which is based the time -series 
development of surface -temperature equations is one expressing the 
transient wall temperature due to an increase of 1° in temperature of 
one wall surface during each interval of time 5. The assumptions 
are made that the Initial temperature is zero and that one wall is 
insulated. From these considerations, F in Duhamel's formula is 
F = t/b and F' = l/S. Substituting this derivative and T from equa- 
tion (B8), with Tp = 1, into equation (B9) and performing the integra- 
tion give the following equation for the temperature at any plane due to 
a surface temperature rise of 1° in the interval 5: 


5 *5 k 6 £3 ( 2 a - 1)5 


• (2n-l)2 If Xt) 
b Z 2 ] 


(Bll) 


The heat flow at any point within the wall due to 1° rise of wall 
surface temperature in the time 5 is obtained by multiplying the con- 
ductivity by the temperature gradient. Differentiating equation (Bll) 
with respect to x and multiplying by K yield 


8 

2 


q = “i. 


IK 



-(2n-l) 


2 *2 
T 
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In order to el imi nate k in the coefficient, substitute for k its 
definition k = K/pc. Then IK/kS = pcl/s. The product pcZ is the 
heat capacity of the -wall per unit area per degree and is represented 
by the symbol G. Si order to obtain the heat flow at the heated sur- 
face, where the heat balance is to be made, let x = l. Hence, the 
instantaneous heat flow due to a uniformly increasing surface tempera- 
ture of 1° per unit of time 5 is simplified to 


q. = 



l 


-(2n-l)2 £| 

(2n - l) 2 


(B12) 


In order to find the average flow of heat over the interval S, 
equation (B12) is integrated with respect to time between the limits 
(t - 5) and t. This integration gives the total heat flow through the 
surface during the interval. On dividing by 5, the average rate of 
flow for the interval is obtained. The result is 


q= ± 9lY~ 

* 2 5 ^i 


(2n - l) 2 


-<2n-l)2 ^ 

^ L 

e 

4 \ 8(2n - l) 1 * 
“L 


-(2n-l)2 id i(t-8) 

X B(2n - l) 4 
4 l d 


(B15) 


The summation 2 — is a constant with the value jt 2 /8. In 

n=l (2n - l) 2 

working with time series only integral increments of time (6, 28, 

58, . . . mS) are used. In equation (B15), therefore, t may have any 
value m8, where m is an integer. For convenience, and to systematize 
results, the following identities are defined: 
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Ao = !ttL2: i—y 

* 2 k5 t3 (2n - 1)4 


-(2n-l) 2 £ i* 
, o ja- 4 12 

a, = it j!_y~ s — 

« 2 k8 n=l (2n - 1)^ 


-2(2n-l) 2 -^6 

p . o° ^ i 2 

A 2 = 'uS-ti t . — 

rt 2 kS n=l (2n - 1H 


-m(2n-l) 2 S- -^-5 
o „ » 4 i2 

A = 4ii_ > e 

( 2n - 1)* 


-(m-l)(2n-l) 2 

— -*-5- 

jc 2 k8 n=l (2n - I) 2 * - 


Equation (Blj) then becomes 


^ " j 2 &(*8 + ^ " ^-1 


for the time interval ending at t = m8. Equation (B15) is used in 
forming the heat balance at the heated vail surface in the derivation of 
wall temperature formulas. 
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APPENDIX C 


ANALYTICAL TEMPERATURE DISTRIBUTION 


Equations (■ 4 - 0 ) are the formulas most convenient for determining 
the temperature of the inside surface or of any plane within the surface. 
For the determination of the temperature of planes within the wall, equa- 
tions (36a) are ordinarily evaluated numerically to obtain 0 for a 
given kt/Z^ and x/l. However, it may be necessary to obtain the tem- 
perature distribution through the wall at some instant of time due to an 
arbitrary surface -temperature history in a form that is analytical with 
respect to x and which can be differentiated or integrated with respect 
to x (for example, in the derivation of a general formula for the thermal 
stress distribution or maximum stress ) . This temperature distribution is 
obtained by substituting equations (36a) into equations (39), substituting 
equations (39) into equations ( 40 ), and collecting terms. The results are 


L i,l 


16 z£ T ( -1) 11 " 1 

3 k5 1 E=I (2n - 1)3 


T-l - 12 . i— T- 


cos 


- «* f] 


-(2n-l) 2 £ 


1 - e 


I2| 


KO n _]_ 


(-l) n_1 cos |(2n - l)|- 


(2n - 1)3 

\2 jt 2 k8 


To - T-, + 


-(2n-l) 2 m -2(2n-l)2 *£ W 

( 2 T X - T 2 )e 4 1 - Tj_e 4 ld 


T i,3 " 


_ _ 00 

16 z£y- 

f 

(-l) n-4 cos J^(2n - 

Si 

k5 H=I 

(2n - 1)3 J 


T3 - t 2 + 


-(20-1)2 B -2(n-l)2 

(-T3 + 2 T 2 - TjJe 1 + ^ 2 T X - T^)e 1 


T]_e 


-3(2n-l) 2 Sf *§- 

4 z 2 
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Note that T ± j_, Tj_ -z )) and Tj_ m can be computed independently. 

The use of these equations involves considerable labor, however, because 
all terms in the summation must be summed in unison. 
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TABIB I.- VALUB8 0E M AND S 5 


\a 5 /l 2 









a \ 

0.01 

0.02 

0.05 

0.1 

0.5 

1.0 

2.0 

5.0 


H 

1 

0.09281491 

0.15125147 

0.20752122 

0.29347746 

0.64728217 

0.85683729 

1.02954122 

1.15145422 

2 

-.01594178 

-.02252780 

-.03560872 

-.05052075 

-.22817193 

-.51142943 

-.82682879 

-1.06920825 

3 

-.02682329 

-.03794254 

-.06002358 

-.08722726 

-.30052433 

-.31662004 

-.20127594 

-.08224362 

4 

-.00957530 

-. 0135^153 

-.02166582 

-.03696546 

-.06405219 

-.02634648 

-.00142616 

-.00000036 

5 

-.00553400 

-.00782681 

-.01320028 

-.02628367 

-.02447704 

-.00223431 

-.00001026 

-.00000000 

6 

-.00373390 

-.00528488 

-.00983262 

-.02025561 

-.00712803 

-.00018948 

-.00000007 


7 

8 
9 

-.00274116 

-.00212389 

-.00170888 

-.00589295 

-.00304410 

-.00249269 

-.00809228 

-.00695594 

-.00608368 

-.0157961c 

-.00207577 

-.00060449 

-.00017604 

-.00001607 



-■ 0123^881 

— .IXJWUJJO 
-.00000012 



10 

-.00141415 

-.00211866 

-.00535620 

-.00733254 

-.00005126 

-.00000001 



11 

-.00119613 

-.00185531 

-.00472749 

-.00588550 

-.00001493 




12 

-.00105005 

-.00166301 

-.00417648 

-.00459660 

-.00000435 




13 

-.00090099 

-.00151704 

-.00369096 

-.00559308 

-.00000127 




Ik 

-.00079666 

-.00140173 

-.00326234 

-.00280743 

-.00000057 




15 

-.00071651 

-.00130718 

-.00288362 

-.00219357 

-.00000011 




16 

-.00064769 

-.00122688 

-.00254890 

-.00171393 

-.00000005 




17 

-.00059289 

-.00115664 

-.00225306 

-.00133917 

-.00000001 




18 

-.00054714 

-.00109385 

-.00199156 

-.00104635 





19 

-.00050882 

-.OOIO3656 

-.00176041 

-.00081756 





20 

-.00047652 

-.00098373 

-.00135610 

-.00063880 





21 

-.00044898 

-.00095449 

-.00157549 

-.00049912 





22 

-.00042534 

-.00088826 

-.00121584 

-.00038998 





23 

-.00040474 

-.00084473 

-.00107473 

-.00050471 





24 

-.00038678 

-.00080337 

-.00094999 

-.00023809 





25 

-.00037081 

-.00076453 

-.00083973 

-.00018602 





26 

-.00035655 

-.00072753 

-.00074227 

-.00014335 





27 

-.00034370 

-.00069238 

-.00065612 

-.00011357 





28 

-.00053195 

-.00065894 

-.00057997 

-.00008874 





29 

-.00032118 

-.00062716 

-.00059695 

-.00051263 

-.00006933 





50 

-OOO31123 

-.00045315 

-.00005417 






6 . . ! 

1 ' 

1.00000000 

0.99999990 

0.99956261 

0.96873107 

0.69945338 

0.45623848 

0. 24814437 

0.09999955 

2 

-.0000C&19 

-.00019242 

-.02166158 

-.12553131 

-.48642965 

-.41618822 

-.24630212 

-.09999906 

3 

-.00002493 

-.00553067 

-.06963676 

-.17636952 

-.15098863 

-.03665380 

-.00182900 

— ocoooo45 

4 

-.00033480 

-.01195026 

-.09012900 

-.14882949 

-.04396967 

-.00310842 

-.00001315 

-.00000000 

5 

-.00144190 

-.02167827 

-.08988005 

-.11749067 

-.01280454 

-.00026361 

-.00000009 


6 

-.00344274 

-.02904777 

-.08284990 

-.09193126 

-.00372885 

-.00002236 

-.00000000 


7 

-.00595604 

-.03361672 

-.07435643 

-.07184409 

-.00108589 

-.00000190 



8 

-.00854572 

-.0359962 6 

-.06609621 

-.<5613647 

-.00031622 

-.00000016 



9 

-.01092536 

-.03687745 

-.05854671 

-.04386265 

-.00009209 

-.OOOOOCCd 



10 

-.01296009 

-.03679326 

-.05179174 

-.03427135 

-.00002682 

-.00000000 



11 

-.01461454 

-.03611078 

-.04579381 

-.02677771 

-.00000781 




12 

-.01590635 

-.03307182 

-.04048318 

-.02092260 

-.00000227 




13 

-.01687671 

-.03383231 

-.03578599 

—01634774 

-.00000066 




14 

-.01737368 

-.05249148 

. -.03163302 

-.01277320 

-.00000019 




15 

-.01804424 

-.03111188 

-.02796174 

-.00996026 

-.00000006 




16 

17 

-.01833036 

—01846822 

-.02973243 

-.02837688 

-.02471644 

-.02184780 

—00779602 

-.00609293 

-.00000002 

-.00000000 




18 . 

-.01848810 

-.02705930 

-.01931206 

-.00476067 





19 

-.01841482 

-.02578770 

-.01707064 

-.00371972 





20 

-.01826876 

-.02456606 

-.01508936 

-.00290638 





21 

-.01806620 

' -.02339609 

-.01333804 

-.00227088 





22 

-.01782041 

-.02227782 

-.01178998 

-.00177434 





2? 

-.01754197 

-. 023 . 21,042 

-.01042159 

-.00138637 





24 

-.01723931 

-.02019254 

-.00921205 

-.00108323 





23 

- .01691916 

-.01922244 

-.00814285 

-.00084638 





26 

-.OI658696 

-.01829829 

-.00719776 

-.00066131 





27 

-.01624686 

-.01741809 

-.00636236 

-.00051671 





28 

' - .01590228 

-.01658000 

- .00562392 

-.00040373 





29 

-.01555588 

-.01578204 

-.00497119 

-.00031545 





30 

-.01520973 

-.01502235 

-.00439422 

-. 0002 46 MS 








‘For additional value of. ItB/l 2 between 0.1 and 1.0, Bee table XI. 
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TABLE II— -VALUES CF M AND 6 USED IN ILLUSTRATIVE EXAMPIES 


'\k 5 / 2 2 
m \ 

0.00980 

0.08820 

0.17640 

0.35280 

0.70560 


M 

i 

0.091882 

O.275619 

0.389750 

0.549116 

0.751602 

2 

-.015782 

-.047357 

-.071018 

-.144143 

-.352559 

5 

-.026554 

-.080945 

-.133036 

-.242444 

-.330551 

4 

-.009479 

-.032753 

-.065871 

-.094474 

-.056482 

5 

-.005478 

-.022958 

-.042292 

-•039558 

-.009904 

6 

-.003696 

-.017994 

-.027361 

-.016565 

-.001737 

7 

-.002714 

-.014408 

-.017705 

-.006936 

-.000305 

8 

-.002103 

-.011580 

-.011457 

-.002905 

-.000053 

9 

-.001692 

-.009314 

- . 007414 

-.001216 

-.000009 

10 

-.ooi 4 oo 

-.007493 

-.004798 

- 000509 

-.000002 

11 

-.001184 

-.006027 

-.003105 

-.000213 


12 

-.001019 

-.004848 

-.002009 

-.000089 


13 

-.000891 

-.003900 

-.001300 

-.000037 


14 

-.000789 

-.003137 

-.000841 

-.000016 


15 

-.000707 

-.002524 

-.000544 

-.000007 


16 

-.000640 

-.002030 

-.000352 

-.000003 


17 

-.000585 

-.001633 

-.000228 

-.000001 


18 

-.000539 

-.001314 

-.000148 

-.000000 


19 

-.000501 

-.001057 

-.000095 



20 

-.000469 

-.000850 

-.000062 




8 

1 

0.999992 

O.992899 

0.943645 

0.804779 

0.580382 

2 

.000008 

-.098507 

-.277731 

-.448794 

-.474632 

3 

-.000021 

-.157664 

-.233674 

-.206910 

-.087207 

4 

-.000290 

-.141628 

-.152510 

-.086651 

-.015291 

5 

-.001288 

-.116097 

-.098715 

-.032684 

-.002681 

6 

-.003146 

-.093525 

-.063879 

-.015194 

-.000470 

7 

-.005529 

-.075429 

-.041337 

-.006362 

-.000082 

8 

-.008028 

-.060636 

- . 026749 

-.002664 

-.000015 

9 

-.010355 

-.048778 

-.017309 

-.001116 

-.000002 

10 

-.012371 

-.039239 

-.011201 

-.000467 

-.000000 

11 

-.014028 

-.031565 

-.007248 

-.000196 


12 

-.015340 

-.025391 

-.004690 

-.000082 


13 

-.016338 

-.020426 

-.003035 

-.000034 


14 

-.017068 

-.016431 

-.001964 

-.000014 


15 

-.017574 

-.013217 

-.001271 

-.000006 


16 

-.017897 

-.010632 

-.000822 

-.000003 


17 

-.018070 

-.008553 

-.000532 

-.000001 


18 

-.018124 

-.006880 

-.000344 

-.000000 


19 

-.018081 

-.005535 

-.000223 



20 

-.017966 

-.004452 

-.000144 
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(a) l/2-inch-thlck wall. 

Figure 1.- Example 1. Temperatures of copper-wall surfaces. Adiabatic- 
wall temperature varies linearly from 0° to 10,000° in 10 seconds; 
h = 100 Btu/(hr) (sq. ft)(°F). 
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(b) 5-inch-thick wall. 


Figure 1.- Concluded. 
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Figure 2.- Example 2. Temperatures of l/2-inch copper wall heated 
according to assigned history of h and T aw . 
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Figure 3« - Example 3* Temperatures of 3~inch copper wall heated 
according to assigned history of h and T av . 
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Present method; 
(8=1/8 sec for 2 s6c; 
18 = 1/2 sec thereafter 
(8=1/8 sec for I sec; 
18 = 1/2 sec thereafter 
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Figure 4.- Example 4. Temperatures of l/2-inch copper wall after appliea 
tion of 5,000 °F jump in gas temperature,, h = 100 Btu/(hr)(sq. ft)(°F) 
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Exact theory for infinite conductivity 
8 = I sec 

S = 2 sec ^ present method; equation (13) 
8 = 5 sec 
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Figure 5-- Example 5* Temperatures of l/l6-inch Inconel vail. 
T aw = 5,000 °F; h = 50 Btu/ (hr) (sq_ ft)(°F). 


NA.CA TN 4105 


54 


Heat 

flow, 

Btu/(hr)(sq ft) 



Time, sec 


Figure 6.- Example 6. Rate of heat flow into heated (outer) surface 
of l/2 -inch copper wall computed^ from temperature history of outer 
surface and from temperature history of inner surface. 
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